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The KamLAND and Borexino experiments have detected electron antineutrinos produced in the 
decay chains of natural thorium and uranium (Th and U geoneutrinos) . We analyze the energy spec- 
tra of current geoneutrino data in combination with solar and long-baseline reactor neutrino data, 
with marginalized three-neutrino oscillation parameters. We consider the case with unconstrained 
Th and U event rates in KamLAND and Borexino, as well as cases with fewer degrees of freedom, as 
obtained by successively assuming for both experiments a common Th/U ratio, a common scaling 
of Th -I- U event rates, and a chondritic Th/U value. In combination, KamLAND and Borexino can 
JlJ , reject the null hypothesis (no geoneutrino signal) at 5a. Interesting bounds or indications emerge 

^ ■ on the Th -I- U geoneutrino rates and on the Th/U ratio, in broad agreement with typical Earth 

model expectations. Conversely, the results disfavor the hypothesis of a georeactor in the Earth's 
I core, if its power exceeds a few TW. The interplay of KamLAND and Borexino geoneutrino data is 

highlighted. 
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(D I. CONTEXT 

The study of low-energy electron antineutrinos (Fg) emitted in the decay chains of uranium (U), thorium (Th), 
and potassium (K) in the Earth's interior — the so-called geoneutrinos — is raising increasing interest in both particle 
' physics and Earth sciences, as recently reviewed in [U, 

■ From the viewpoint of particle physics, there has been dramatic progress in understanding the flavor evolution and 
^ I oscillations of neutrinos [3], and in refining their low-energy, low-background detection techniques. In particular, the 
. detection of the global Th -I- U geoneutrino flux in the KamLAND [J [H] and Borexino @ experiments through the 
inverse beta decay reaction 



-p^n + e+ iE^> 1.806 MeY) , (1) 



. . , represents a milestone in this field. ^ In perspective, measurements at different locations might constrain the relative 
^ ' Th and U abundances in different reservoirs (e.g., crust versus mantle), especially if some directional sensitivity could 
k>( [ be achieved 

5_j ■ From the viewpoint of Earth sciences, the heat-producing elements U, Th and K, despite their relatively low natural 
d ' abundances, bear on outstanding and debated issues. Their global amounts should reflect different condensation 
histories in the primitive solar nebula, which presumably led to partial escape of the moderately "volatile" K and 
to complete capture of the "refractory" U and Th. Therefore, the mass abundance ratio of Th and U in the Earth 
(Th/U hereafter) is expected to be the same as in the most pristine meteorite samples, the so-called carbonaceous 
chondrites 

Th 

-jj- ~ 3.9 (chondritic estimate) . (2) 

The radial distribution of U, Th and K should instead reflect subsequent Earth differentiation processes, as these 
elements are both "lithophile" (preferring mantle and crust silicates to core metals) and "incompatible" (preferring 
crust melts to mantle residues). Within the mantle, slow convection processes may have further redistributed the 
radiogenic elements in several possible ways Q. 



^ Geo-u from K decay are below threshold for the reaction JTJ. 
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In this field, connecting geophysical quantities to particle physics observables requires some modeling of the U, 
Th, and K distributions, on both planetary and local scales. For instance, a correlation is expected between the 
radiogenic heat production rate H and the geoneutrino event rate R from Th + U sources, within large uncertainties 
induced, e.g., by the unknown distribution of radiogenic elements in the mantle. Various {H, R) correlation plots have 
been discussed in 1] under rather general and plausible assumptions [including the estimate in ([2])] and conservative 
uncertainties. In particular, the results of [ijj for KamLAND (see Figs. 23 and 30 therein) can be approximated as 

^^^^ ^ ~ (1.11 ± 0.14) X ^^^^ ^ - 25.0 (KamLAND, 6 TW < H < 40 TW) , (3) 
TW ^ ' TNU ^ - - ; > \ J 

where 1 TNU (terrestrial neutrino unit) corresponds to lO'^^ events per target proton per year, and the quoted 
error provides a sort of "maximum allowed range" for the (i/, R) correlation band. The above estimate holds for 

> 6 TW, which provides a total KamLAND signal in excess of the "guaranteed" minimum contribution from the 
Earth's crust: i?(Th + U) > 24 TNU [1]. On the other hand, Eq. © is not apphcable beyond the "fully radiogenic" 
limit if(Th + U) ~ 40 TW [l|, which, adding an estimated potassium contribution H{K) ~ 5 TW f^, would saturate 
the global Earth's heat flow, ff® - 45 TW [IQ. 

Positive correlations are also expected among the expected event rates in different experiments, since they probe 
the same geoneutrino sources, although weighted differently by the inverse square law for the fluxes. In particular, it 
is rather plausible to assume that KamLAND (KL) and Borexino (BX) probe the same average Th/U ratio, so that 



i?(U)BX i?(U)KL 



i?(Th)BX i?(Th)KL 



(4) 



Even if the primordial proportions of Th and U in the Earth were different from the estimate in ([2]), the known 
geochemical similarity of Th and U in different reservoirs would support the above assumption. 

Correlations among absolute rates (rather than ratio of rates) may be more model-dependent. For instance, a 
comparison of the KamLAND (KL) and Borexino (BX) uranium event rates estimated in [l| for a wide range of 
admissible Earth models (see Table 12 and Fig. 24 therein) suggests an approximate scaling law, 

i?(U)BX^1.15i?(U)KL , (5) 

within a relatively small spread (< 5%) in the coefficient. A scaling coefficient > 1 is to be expected, since Borexino 
probes a thicker crust than KamLAND. Its spread, however, may actually be larger than 5%, since the local (and 
sizable) crust contributions at the two sites are not expected to have significant covariances, thus reducing the overall 
correlation [To| . Nevertheless, within the large uncertainties affecting current geoneutrino event rates, an approximate 
scaling assumption as in Eq. ([S]) can still provide a useful guidance in the data analysis, and will be used later. 

Summarizing, the most general analysis of the available KamLAND and Borexino data involves four geoneutrino 
degrees of freedom {Nu = 4), namely, the thorium and uranium event rates in the two experiments. This parameter 
space can be reduced by assuming either of the two Eqs. Il4|5p . If the estimate in ([2]) is also imposed, a single degree 
of freedom remains {Nd ~ 1). The main purpose of this work is to analyze geoneutrino data, and to discuss their 
implications, in cases with Njj = 4, 3, 2, and 1. A special case is also discussed, where an additional degree of freedom 
is provided by the unknown power of a hypothetical georeactor in the Earth's core [Tll |. 

Our paper is structured as follows. In Section II we discuss our approach to the geoneutrino analysis, with particular 
attention to the theoretical and experimental energy spectra and their uncertainties. Since low-energy neutrino 
oscillations — constrained by solar and long-baseline reactor data — affect the extraction of geoneutrino signals, the 
marginalization of the three-neutrino oscillation parameters in a global fit is also discussed. In Section III we show 
the results of our analysis in four relevant cases [Nd = 4, 3, 2, and 1), and discuss their implications. In general, 
we find results in agreement with typical Earth model expectations, although often within large uncertainties. An 
interesting interplay between KamLAND and Borexino data emerges in all cases with Nd < 3. Finally, we show that 
the georeactor hypothesis is disfavored by the data. Our conclusions are summarized in Section IV. 



II. METHOD 



In this Section we describe some aspects of our analysis, concerning the geoneutrino and reactor energy spectra, the 
marginalization of oscillation parameters, the input data from KamLAND and Borexino, and the relevant geoneutrino 
degrees of freedom under increasingly restrictive assumptions about the relative event rates. 
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FIG. 1: Geoneutrino spectra computed in this work, as a function of the neutrino energy E^. Upper panel: spectra from U 
and Th decay chains (solid and dashed curves, respectively). Middle panel: spectra multiplied by the cross section for inverse 
beta decay. Lower panel: observable event spectra, including typical energy resolution effects. The results refer to a Fe flux of 
10^/cm'^/s (all panels), and to 10"^^ target protons (lower panel). 



A. Geoneutrino energy spectra 

The reaction ([T]) allows to detect (a fraction of) the 4 X/g and the 6 produced in the decay chain of the ^^^Th and 
the 238u nuclei, respectively, as described at length in 1], to which we refer the reader for details. We perform our 
own calculation of the corresponding energy spectra, based on experimental nuclear data from fT3| and theoretical 
inputs from [l3| . 

Our results are reported in Fig. 1. The upper panel shows the computed energy spectra X{E) from the Th and U 
decay chains, assuming in both cases a reference 17^ flux of 10^/cm^/s. The spectra are normalized to unit area for 
>0 (lower range not shown). From our spectra we estimate that the flux fraction above the 1.806 MeV threshold is 
0.151/4 (Th) and 0.384/6 (U). The middle panel of Fig. 1 shows the spectra times the inverse beta decay cross section 
[3]. We estimate the average cross sections as crxh = 12.9 x 10~^^ cm^ and au = 40.5 x 10~*^ cm^, corresponding to 
the event rates i?(Th) = 4.07 TNU and i?(U) = 12.8 TNU for a reference Ue flux of 10*^/cmVs. Since 1 kg of natural 
Th emits 16.2 x 10^ Ve/s from ^^^Th decays, and 1 kg of natural U emits 74.1 x 10^ Ve/s from ^^®U decays [l|, the 
event rates and the natural mass abundances of Th and U for a given source are related by 

The lower panel of Fig. 1 shows the typical effect of finite energy resolution in liquid-scintillator detectors such as 
KamLAND and Borexino. The event rate spectra appear to be significantly smeared out, in the range E,y ~ 1.7- 
3.5 MeV. Separation of Th and U contribution requires, in principle, an accurate determination of relative event rates 
above and below E^ ~ 2.5 MeV. Experimentally, one does not measure E^, but the observable "prompt" energy Ep 
associated to the final-state positron and its annihilation 



Ep = + 2me ^ E^- (m„ - rrip) + nie :^ E^ ~ 0.782 MeV . (7) 

Accounting for smearing effects on Ep, approximate geoneutrino energy windows for Th and U events are then 
Ep{Th) e [0.9, 1.7] MeV and Ep{\J) e [0.9, 2.6] MeV, respectively. 
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B. Reactor energy spectra 

Concerning KamLAND, we have described our approach to the reactor spectra and data analysis in previous works 
0, [iB4i3j to which we refer the reader for details. Concerning Borexino, we integrate available information on 
position, type and average power of European and world reactors from several public sources, see e.g. jl^, and adopt 
the typical power fractions of fuel components as suggested in Q. Since the oscillated reactor spectra at Borexino 
are largely averaged over time and over many (and very long) baselines — besides being smeared by energy resolution 
effects — more accurate information is not really needed for our purposes. 



C. Marginalization of oscillation parameters 

In KamLAND and Borexino, geoneutrino events can be distinguished by reactor and background events only 
statistically, on the basis of their different energy spectra. The reactor spectra are significantly affected bythree- 
neutrino oscillations governed by the squared mass gap Srn^ — — ni\ and by the mixing angles 6*12 and ^13 [3|. The 
pattern of 5m^-driven oscillations clearly emerges in KamLAND reactor spectra after pathlengths L ^ O(IO^) km @, 
while the pattern is largely averaged out in the case of Borexino, where L ~ 0(10'^) km [6j. Complete averaging of 
oscillations can be assumed, to a good approximation, for geoneutrinos [ij . 

The statistical separation of geoneutrino, reactor and background spectra in KamLAND and Borexino depends 
thus on the oscillation parameters {5rn? , 0i2, ^13) which, in turn, are also constrained by solar neutrino data Q. We 
perform a combined analysis of KamLAND and Borexino data, together with all solar neutrino data, updating our 
previous work In particular, we include the latest Gallium experiment event rates [13, H!! and the low-energy 

threshold data from the Sudbury Neutrino Observatory (SNO) [22|. We also consider both low- and high-metallicity 
options for the Standard Solar Model [1^. The statistical analysis involves a 7-dimensional manifold, spanned by 

{(5m2, 012, 013; i?(Th)KL, i?(U)KL, i?(Th)BX, i?(U)Bx}, (8) 

plus a number of nuisance parameters which account for systematic uncertainties via the pull method (2^ (see also 
the next subsection). 

In the global data fit, the marginalization of the four geoneutrino rates R provides bounds on the low-energy oscilla- 
tion parameters (5m? , 0i2, 0i3), which will be discussed elsewhere. In this work we are interested in the complementary 
case, where the marginalization of {5m? , 0i2, ^13) provides constraints on the Th and U event rates in KamLAND and 
Borexino. 



D. KamLAND and Borexino data 



Experimental geoneutrino spectra are usually presented in terms of number of events Ni observed in bins of prompt 
energy (or a related variable) . The total number of events N corresponds to an event rate R at the detector via 

N = eTR, (9) 

where e is the detection efficiency, and T is target exposure in units of 10^^ protons x year, with R expressed in TNU. 

We use the latest KamLAND spectral data after an exposure T — 2.44 x 10'^^ protonsxyear with small-width 
bins in the geoneutrino range Ep < 2.6 MeV [1^. We include the energy-dependent function e{Ep) from Q, which 
implies average detection efficiencies exh = 0.69 and ey = 0.78 for Th and U events, respectively A Poisson 

function is constructed [l^ to account for statistical fluctuations. We introduce five systematic pulls, one for the 
energy scale uncertainty, and four for the normalization of: (i) reactor events; (ii) all events; (Hi) ground-state and 
(iv) excited-state contributions to the ^'^C(a,n)^^0 background. 

Figure 2 shows the experimental KamLAND spectrum in terms of prompt energy, as well as the partial and total 
contributions to the theoretical spectrum (background, reactor signal, and geoneutrino signal) at the best fit in the 
parameter space (|8]). The geoneutrino signal is more pronounced at low energies, consistently with a relatively large 
contribution from Th decay {Ep < 1.7 MeV). Our geoneutrino fit results are very similar to the official ones [5|; 
moreover, we reproduce the fit of [25| in terms of total Th -I- U rate versus Th — U rate asymmetry (not shown). We 
find that the null hypothesis (no geoneutrino signal in KamLAND) is rejected at 2.9cr. 

Concerning Borexino, we use the spectral data for an effective exposure eT = 0.152 x lO^^protonsxyear, and in 
terms of the light yield Y for positron events — approximately equal to Y ~ 500 x Ep/MeV 0. The data from the 
americium-bery Ilium source calibration in [6*1 allow us to improve this approximation for Y , as well as to infer the 
energy resolution width in Borexino. Our Poisson statistics includes two systematic pulls for the normalization of 
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FIG. 2: KamLAND event spectrum as a function of the observed prompt energy Ep. Data points and backround (BG) estimates 
are taken from [1, jl^; note the narrower binning in the geoneutrino energy range Ep < 2.6 MeV [2^. The histogram represents 
our best-fit spectrum, with cumulative contributions from background, plus reactor, plus geoneutrino events. 

all events and of reactor-only events. Since the background is very small in Borexino, its systematic uncertainties are 
negligibile for our purposes. 

Figure 3 shows the experimental Borexino spectrum in terms of light yield, as well as the separate contributions 
to our best-fit spectrum, in analogy with Fig. 2. Note the clear geoneutrino signal, which covers the whole expected 
range Ep € [0.9, 2.6] MeV {Y S [450, 1300]), consistent with a leading contribution from U decay. Also in this case, 
we are able to reproduce quite well the official geoneutrino fit results and plots of ^] (not shown). We find that the 
null hypothesis in Borexino is rejected at 4.1cr. 
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FIG. 3: Borexino event spectrum as a function of the light yield for positron events (p.e.). Data points and background (BG) 
estimates are taken from @]. The histogram represents our best-fit spectrum, with cumulative contributions from background, 
plus reactor, plus geoneutrino events. 
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TABLE I: Summary of adopted degrees of freedom and constraints. 



Nd 


Constraints 




i?(Th + U)KL 


(Th/U)KL 


i?(Th + U)bx 


(Th/U)BX 


4 


None 




free 


free 


free 


free 


3 


(Th/U)BX = 


(Th/U)KL 


free 


free 


free 




2 


(Th/U)BX = 


(Th/U)KL and Rbx = 1.15 i?KL 


free 


free 






1 


(Th/U)BX = 


(T1i/U)kl = 3.9 and Rbx = 1.15 i?KL 


free 









E. Analyses with 4, 3, 2, and 1 degrees of freedom 

After marginalization of the oscillation parameters in (jSj, the most general geoneutrino parameter space is spanned 
by four event rates. We find it useful to define four equivalent geoneutrino degrees of freedom {Nd = 4), namely, the 
total rate and the Th/U ratio probed by KamLAND and Borexino, 

i?(Th + U)KL = i?(Th)KL + i?(U)KL , (10) 

(Th/U)KL = [i?(Th)KL/i?(U)KL] /6.96 X 10-2 , (11) 
i?(Th + U)BX = i?(Th)Bx + i?(U)BX , (12) 
(Th/U)Bx = [i?(Th)Bx/i?(U)Bx] /6.96 x IQ-^ . (13) 

One can reduce the degrees of freedom to Njj = 3 by assuming that KL and BX probe essentially the same Th/U 
ratio, 

(Th/U)Bx = (Th/U)KL . (14) 

If, in addition, the scaling law in Eq. ([5]) is assumed, then A^^i — 2. Finally, if the chondritic Th/U estimate in ([2]) is 
also assumed, then No — 1. These four options, involving an increasing Earth model dependence for decreasing Nd, 
are summarized in Table HI 

A final remark is in order. As discussed in the next Section, the allowed ranges for the KamLAND and Borexino 
geoneutrino degrees of freedom may extend beyond plausible expectations, where the constraints in Eqs. ([5]) and 
(fT4|l are not really justified by any Earth model. Therefore, while the analysis for Nd = 4 is completely general, the 
results of constrained analyses {Nd < 3) must be taken with a grain of salt. 



F. Analysis with an additional degree of freedom: the georeactor 

It has been proposed [ll| that there could be enough uranium in the Earth's core to naturally start a nuclear fission 
chain over geological timescales, with a typical power (at the current epoch) of Pgeo — 3-10 TW [1^. This hypothesis 
is disfavored by various geochemical and geophysical arguments . Particle physics offers an independent probe of 
the hypothesis, since a georeactor would alter the observable energy (and time) spectra of T^e events |27l]. In particular, 
we reported in an analysis of earlier KamLAND data in the energy and time domain, providing an upper bound 
Pgco ^ 13 TW at 95% C.L. Further KamLAND and Borexino data have reduced the upper bound to < 6.2 TW at 
90% C.L. [| and < 3 TW at 95% C.L. il respectively. 

Here we update our previous analysis jl6| , by assuming a contribution from a georeactor at the center of the Earth 
(with unconstrained Pgoo) in the KamLAND and Borexino energy spectra, for each of the four cases in Table HI With 
respect to [loj. the current analysis does not include the event time information, which has not been released by the 
experiments^, 0| . 



III. RESULTS 



In this Section we describe the results of our analysis, in terms of both joint and separate bounds on the Th/U and 
i?(Th + U) variables. 
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FIG. 4: KamLAND (KL) and Borexino (BX) geoneutrino analysis in tiie plane charted by the total rate _R(Th + U) and by 
the mass abundance ratio Th/U. The curves represents la contours (Ax'^ = 1) around the best-fit points (thick dots). From 
top to bottom, the degrees of fi-eedom decrease fi-om Nd = 4 to Nd = 1, as reported in Table ID 



A. Joint la regions for Th/U and R{Th + U) 

Figure 4 shows the la contours (Ax^ 1) in the plane charted by the total event rate i?(Th + U) and by the mass 
abundance ratio Th/U for KamLAND and Borexino. The degrees of freedom decrease from Nd = 4 to Nd — 1 from 
top to bottom, according to the constraints in Table HI 

The upper panel {Nd — 4) shows that both KamLAND and Borexino place upper and lower bounds to the total 
event rate i?(Th + U). These bounds are consistent with typical Earth model expectations, which span the la range 
29-41 TNU for KamLAND and 34-48 TNU for Borexino (see Table 11 in [1|]). However, neither KamLAND nor 
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TABLE II: Best fits and la ranges from the data analysis with degrees of freedom Nd < 4. Event rates R are expressed in 
TNU. Derived or fixed numbers are given in brackets. 



Nd 


R{Th + U)kl 


(Th/U)KL 


i?(Th + U)Bx 


(Th/U)Bx 


4 




oc q+oo 

^o.»_22.g 


66.91^^:^ 


r, 7+20.2 
^•'-2.7 


3 


41.3tl« 


q 1+23.5 
».r_7.4 




[9.11?:^-] 


2 






[5i.7t}^2:^] 


[9.61?-] 


1 


^' • ' -11.2 


[3.9] 


[54.9il||] 


[3.9] 



Borexino can currently determine Th and U separately. In particular, KamLAND is compatible with all events 
being from Th decay (Th/U= oo), while Borexino is compatible with all events being from U decay (Th/U= 0), 
as anticipated in the context of Figs. 2 and 3, respectively. However, a broad range of Th/U values appears to be 
compatible with both KamLAND and Borexino results at la; this range excludes the extreme cases with null Th or 
U signal, and includes the chondritic value Th/U= 3.9. Therefore, it makes sense to reduce the degrees of freedom 
by imposing that the two experiments probe the same Th/U ratio as in Eq. (fT4)) . 

The second panel {Nd = 3) shows the results of such exercise, providing both upper and lower la limits on the 
Th/U ratio, with a best fit which is only a factor of ^^2 higher than the chondritic value. The total rate estimates 
are not significantly altered with respect to the case with Nd — 4. Therefore, under the rather general assumption 
in Eq. (|T4)) . the combination of KamLAND and Borexino data starts to be sensitive to the global Th/U ratio of the 
Earth, although only at the ^Ict level; as discussed below, current Th/U constraints disappear at ~1.5ct. 

The results in the third panel {Nd ~ 2) include, in addition, the approximate scaling assumption in Eq. ([5]). In 
this case, the KamLAND parameters i?(Th + U) and Th/U are conventionally taken as free, while the corresponding 
Borexino parameters are derived (hence the "dotted" BX contour in the panel). In this case, the best fits for the total 
rates are located slightly above the quoted Earth model expectations (29-41 TNU for KamLAND and 34-48 TNU 
for Borexino [l|), with la uncertainties at the level of ^ 30%, dominated by KamLAND data. Concerning Th/U, the 
best fits and la ranges are not significantly altered with respect to the previous case with Nd = 3. 

The comparison of the three panels with Nd = 4, 3 and 2 shows that the current constraints on the total rates and 
on the Th/U ratio are approximately independent: they do not significantly affect each other within present data. 
The weak negative correlation in the la contours reflects the fact the overall rate i?(Th + U) increases somewhat faster 
for larger U contribution as compared to Th contribution, the latter being confined at low energy (see the spectra in 
Fig. 1). 

Finally, the results in the lower panel of Fig. 4 {Nd = 1) include, in addition to the previous constraints, the 
chondritic estimate in ©.^ In this case, the estimated KL total rate is i?(Th + U)KL = 47. 7± 11.2 TNU, the BX total 
rate being a factor '--^I.IS higher by construction, i?(Th + U)bx = 54.9 ± 12.9 TNU. These results show a preference 
for Earth models with relatively high expectations in Th and U contents, although within very large uncertainties at 
present. We report in Table HIl a numerical summary of the la ranges for the total rate and Th/U ratio, in each of 
the four cases considered. 

In the most constrained case {Nd — 1), where the total rate error is reduced to ^ 23%, it makes sense to infer 
indications about the associated radiogenic heat _ff(Th + U) via the approximate {H, R) correlation in Eq. ([3]). By 
fixing the KL rate at its central value, i?(Th+U)KL — 47.7 TNU, one would obtain an allowed range iJ(Th+U) ~ 21- 
35 TW, somewhat above the plausible expectations of 14-18 TW (see Fig. 23 in However, including the la rate 
uncertainties, the allowed range in significantly enlarged: (Th + U) ~ 10-49 TW. The upper value is not particularly 
meaningful, being larger than the "fully radiogenic" limit of ~40 TW. The lower value of ^ 10 TW, however, exceeds 
the "guaranteed" contribution from Th and U in the crust (~6 TW pj), and suggests, indirectly, the presence of an 
additional contribution from a different reservoir — which can be naturally identified with the mantle. 

In conclusion, the combination of KamLAND and Borexino data brings to surface some intriguing — although still 
statistically weak — pieces of information: {i) preferred Th/U values in broad agreement with chondritic expectations; 
(m) slight preference for relatively high Th and U contents in the Earth; and {Hi) hints of a mantle contribution 
to the total geoneutrino signal. We remark that these indications emerge only at the ^ la level from the current, 
low-statistics data samples. 



^ In the panel, the Th/U coordinates of KL and BX are slightly displaced from 3.9 for the sake of clarity. 
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FIG. 5: KamLAND and Borexino analysis: constraints on the Th/U ratio in terms of standard deviations, A^,,- = ^ Ax^ 



B. Separate bounds on Th/U and Jl(Th + U) 

In this section we discuss the separate projections of the previous results onto the variables Th/U and i?(Th + U), 

in terms of standard deviations from their best fit {Na — yAx^)- 

Figure 5 shows the constraints on the Th/U ratio from KamLAND and Borexino. In the upper panel {Nu = 4), 
one can see at a glance that KamLAND and Borexino data place lower and upper \a limits on Th/U, respectively, 
but have no significant sensitivity in the opposite directions. The bounds are still statistically weak, as they vanish 
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FIG. 6: KamLAND and Borexino analysis: constraints on the total rate -R(Tli + U) in terms of standard deviations A^,,- 



at '^1.3(7 in KamLAND and at '^l.Gcr in Borexino. The middle panels show that KamLAND and Borexino provide 
joint limits on Th/U at the Icr level, vifith no significant variation between the cases with Nu = 3 and = 2. The 
lower panel is just a "Dirac delta" at Th/U= 3.9. 

Figure 6 shows the constraints on the total event rate i?(Th+U) in KamLAND and Borexino. The two upper panels 
shows that the null hypothesis of no geoneutrino signal is rejected at 2.9(t in KamLAND and at 4. Icr in Borexino, 
almost independently of the Th/U constraint in Eq. (jl4l) . The two lower panels show the significant error reduction 
induced by the scaling law assumption in Eq. ([S]), which strengthens the null hypothesis rejection at the 5a level. By 
comparing the two lower panels, it appears that the chondritic assumption in ^ has a minor impact on combined 
rate constraints. 
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In conclusion, current indications on the Th/U ratio emerge only at the level of ^ la, while global geoneutrino 
signals emerge at > 3a, reaching ^5a in constrained combinations. Bounds on Th+U rates and on the Th/U ratio 
are largely uncorrelated. Detailed la ranges for both free and derived parameters are reported in Table HIl 

C. Bounds on the power of a hypothetical georeactor 

We have repeated the various analyses by including a spectral contribution from a hypothetical georeactor with 
unknown power Pgoo, located at the center of the Earth. The results disfavor this hypothesis in all cases, if Pgco 
exceeds a few TW. In the most general case (unconstrained geoneutrino and georeactor event rates), we find that 
KamLAND and Borexino place the 2a upper limits Pgeo ^4.1 TW and < 6.7 TW, respectively. In combination, the 
joint KamLAND+Borexino limit reads 

Pgco < 3.9 TW at 2a (< 5.2 TW at 3a) , (15) 

almost independently of the chosen geoneutrino degrees of freedom. This is to be expected, since the georeactor 
spectrum extends well above the geoneutrino energy range. 

Our combined constraints on Pgeo appear to be dominated by Borexino data, but do not improve upon the official 
Borexino limit (quoted as Pgoo ^ 3 TW at 2a [^), presumably because we cannot include information in the time 
domain, which is currently unpublished. Our analysis of older KamLAND data has indeed shown that the time 
spectra can add significant constraints to Pgoo. 

IV. CONCLUSIONS 

We have performed a detailed analysis of current geoneutrino events from Th and U decay chains as detected 
in KamLAND [5|] and Borexino |6i|, within the more general context of low-energy neutrino oscillation data from 
long-baseline reactor and solar sources, and of a broad range of Earth model expectations taken from Ij. 

The relevant parameter space is spanned by the total Th+U event rates and the Th/U ratio in the two experiments, 
while the oscillation parameters {dm^, 612, ^13) are marginalized away. Constrained analyses with fewer degrees of 
freedom are obtained by successively assuming for both experiments a common Th/U ratio [Eq. (jl4p ]. a common 
scaling of Th + U event rates [Eq. ([5|)], and a fixed (chondritic) Th/U value [Eq. ([2])]. Cases with a hypothetical 
georeactor, involving an additional degree of freedom, are also considered. 

The results are in agreement with typical Earth model expectations, although within still large uncertainties. The 
global Th+U geoneutrino signal emerges at 2.9(t and 4.1cr in KamLAND and Borexino, respectively, and can reach the 
overall 5(j level in combination. The data disfavor the hypothesis of a georeactor, and limit its power to Pgeo ^ 3.9 TW 
at 2a (or Pgoo < 5.2 TW at 3a). 

Weaker — but potentially interesting — pieces of information emerge at the ~ la level, including: {i) preferred Th/U 
values in broad agreement with chondritic expectations; (m) slight preference for relatively high Th and U contents 
in the Earth; and [iii) possible hints of a (mantle) contribution in excess of the "guaranteed" signal from the crust. 
Significantly higher statistics, possibly from new large- volume detectors and with some directional sensitivity, will be 
needed to promote these intriguing indications to more robust signals improving our understanding of the Earth's 
interior. 
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